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1, INTRODUCTION 

Through of stages of delivery electrical power, which are the generation, transmission and 
distribution systems using overhead lines the isolation system exposed to many stresses electrical, 
mechanical and environmental stress. A reliability of power system might be disrupted by reinforcement of 
electric field along of the insulators [1]. Damage to the outdoor electrical insulators resulting of the electrical 
fields could be increased especially near to electrodes exponentially with time [2], [3]. 

The Silicone rubber insulators have some of the characteristics better than some insulators such as 
hydrophobicity, possess higher impact strength, and compact design, so their application on power 
transmission system tends to increase rapidly. Electric fields distribution, observed on insulators used in 
overhead lines, are the most important issue for distribution systems [2]. Field stress is a complex issue due 
to diverse reason such as different contaminations, a shape of insulators, and different stresses on the 
insulator and unknown impact of moisture on the contaminations [4]. The efficiency of insulators under 
different stresses could be tested by analysis of electrical field on their surface. On the other hand, the field 
stress on insulators can cause power failure at a very far distance. The field stress on insulators is still a 
damage threat to the safe process of a power distribution system, which might led to occur of flashover [3]. 
It is generally believed that electric field is becoming ever more critical to investigate the performance of 
high voltage insulators [5]. 

In this article, the electric field using 3D finite element method (FEM) using COMSOL 
Multiphysics 2.5 software based on replacement and modify of insulators profile along a standard 11 kV 
silicone rubber insulator has been computed. This studying took to account inclination of angle 8, shed radius 
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R, shed length LSH, and distance from shed to ground terminal L on the 11kV silicon rubber insulator 
Moreover, a suitable solution to reduce of electric field were suggested. This can be considered as an 
approach to enhance & optimized the insulation performance, similar to other dielectric materials i.e. gases 
and liquid, where researchers try to find best possible candidate as for an alternative [6], [7]. 


2. MATHEMATICAL CHARACTERIZATION 
2.1. Calculate of E-Field and Potential Voltage Distribution 

One easy method for field stress determination on insulators surfaces is determine voltage 
distribution [4]. Then, the E-field distribution is calculated by minus gradient of voltage distribution. 
In electrostatic field case, E-field distribution can be defined as follows [5]: 


E=—VV (1) 
Based on Maxwell’s equation, 

VE =p/e (2) 
Where E is the electric field, V Apply voltage, p Resistivity Q/m and ¢ Material permittivity 

According to relative permittivity (er), the isotropic material was classified. Relative permittivity or 
dielectric constant (er) linked the material permittivity (€) with the permittivity of vacuum (0) as the 
following equation: 

f=" 2,4, ='3.85 * 10 = (3) 
For obtaining Poisson Equation place (1) to (2) 

eV7V = —p¢ (4) 
As shows in (4) becomes Laplace equation when p = 0 as follows: 

eV*V =0 (5) 
2.2. Analysis of Finite Element Method (FEM) 

Many software like ANSYS, ADINA, LS-DYNA, and COSMOL Multiphysics 5.2 use FEM to 
compute of structures properties, three stages in each FEM processing are pre-processing and manipulation 
of data, compute of FEM case, and post-processing of final results. FEM has been vastly used in electric, 
magnetic, and thermal application as well as hydraulic applications since the late1970s [8]. 


Assuming that the area under research does not include charges for any space and surface. The 2D 
functional F(u) in the coordinates Cartesian system as in Figure | can be expressed as follows [9]: 


F(u) = =f, le, (“y + &y (“) | dxdy (6) 


Where €, 1s X- Components of material permittivity, €, 1s Y-Components of material permittivity in the and 
u = Potential voltage. 
In the isotropic permittivity case (€ = €, = €,) and 


F(u) = =f, € (4) + (*) | dxdy (7) 


If influence dielectric loss on the distribution of electric field took to account, the complex function 
F(u) should be written as 


F(u)* = =f, WE (€ — jetgd) (2) + (=) ] dxdy (8) 
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Where @ is angular frequency, tgo = the tangent of the insulator loss angle and u* is the complex potential 
voltage. Assuming there 1s linear diversity of the voltage in each sub-domain Se. 


Ue(X, VY) = Bey + AenX + Aozy ; (e=1,2, 3,....ne) (9) 


Where is @,, 18 coefficients of computational for element el, @,2 coefficients of computational for element 
e2, Ao; coefficients of computational for element e3, u,(x, y) potential voltage of any arbitrary dot inside 
each sub-domain Se and Se is number of elements of triangle. The electric potential at the node in the 
network can be calculated by the following function F(u), that 1s: 


OF (uj) 
Ou; 


=0;i=1,2,....np (10) 


Where np is stands for the total number of nodes in the network. 


3. CHARACTERIZATION FOR MODEL 
3.1. Properties and Materials of the Model Insulator 

Four sheds 11kV polymeric (silicone rubber SIR) insulator was selected to compute of the electric 
field and potential voltage distribution. The polymeric insulator made up of three main materials as illustrated 
in Figure. 1:(1) Loadbearing and insulating material a fibre reinforced plastic (FRP) as the rod, (11) Insulating 
material silicone rubber as weather sheds that cover the to protect it from the vandalism(Polymeric Housing), 
(111) Loadbearing material metal end fittings that made of forged steel to withstand the mechanical 
stress [10], [11]. The dimensions of insulator that showed in Figure | are presented in Table 1. 








End fittings 


Weather 
sheds 


FRP rod 


Figure 1. An 11kV SIR insulator under study:(a) A 3D insulator model, and (b) practical insulator 


Table 1. Details the Dimensions of Insulator 


Parameter Length in (mm) 
Creepage distance 375 
End fitting separation 175 
Core diameter 18 
Diameter of trunk 28 
Diameter of shed 90 
Separation of shed 46 


The properties of material selected for modelling of insulator in this research are demonstrated in 
Table 2 [10]. 
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Table 2 Properties of Material for Modelling of Insulator 
Conductivity, o(uS/m) 


Material Relative Permittivity, ¢, 
FRP core 7.1 0.1 x 1078 
SIR 4.3 0.1 x 10-8 
Forged steel 1 5.9 x 1013 
Layer of pollution 80 0.6 
Air 1 0.1 x 10~? 


3.2. Geometry and Physics of Insulator 
An eleven kV at 50 Hz is energized on the upper terminal while the lower terminal connected 


directly to the ground (0 kV). The material properties as was mentioned before having applied to the domains 
of the Three-dimensions of the silicone rubber insulator for FEM simulation using COMSOL Multiphysics 


5.2 software as illustrated in Figure 2(a) [13], [14]. To boost simulation accuracy the mesh element along the 
leakage distance of the insulator will refine manually [12], [15]. The insulator and its surrounding air medium 


represented in Figure 2b. 
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Figure 2. (a) A 3D SIR insulator model in pre-simulation (b) Meshed model 


4. SIMULATION OF THE MODEL 
The electrostatic formulation of AC/DC physics was used for all simulations. Effect of variation of 


shed radius at R= 5, 10, and 15 mm), shed’s inclination angle at 8 = 5, 10, and 15 mm, shed’s length at Lsh = 
45, 65, and 85 mm, diameter of core at r=28, 33,and 38, and distance from the shed to the end fitting at H= 
20, 30, and 40 mm proposed on the electric field and electric potential distribution has studied in detail as 


showed in Figure 3. 





IU 


(a) Change of the shed radius 
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(e) Change of fitting distance 


Figure 3. The proposed profiles for SiR insulator in clean and pollution conditions 


5. RESULT OF SIMULATION 

Potential voltage and electrical field distribution in the 11-kV clean polymeric insulator in the 
normal case is presented in Figure 5. In the condition of the clean insulator, a small current called as 
displacement current is flowing over the surface. According to the self-capacitance, the voltage for a clean 
insulator surface is distributed. 


5.1. Effect of Change of the Shed’s Radius 

Electric Potential and electric field distribution along silicon rubber insulator is shown in Figure 
4(b). The shed’s radius on the insulator was varied to simulate and study its effect in electric potential 
and field. 


Indonesian J Elec Eng & Comp Sci, Vol. 14, No. 2, May 2019 : 618 — 627 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 623 


A slight change in the electric field and potential in the surround of ground electrode region of the 
dielectric while in the surround of high voltage area the field intensity and voltage increase with the increase 


of the shed’s radius as demonstrated in Figure 5. 
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Figure 4. (a) Electric potential and (b) Electric field distribution under clean condition 
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Figure 5. (a) Electric potential and (b) Electric field with different shed’s radius 


5.2. Effect of Change of the Shed’s Inclination Angle Diameter 


The shed's inclination angle diameter of the 11kV RIS insulator under study is 5mm as shown in 
above Figure 3(b). The parameter is varied with values of 5, 10 and 15 to study the effects of the electric field 
stress. The electric potential and electric field stress along the surface of the insulator is shown in Figure 6. 

The increase of electric field stress and the voltage near to the electrodes with the increase of the 
shed’s inclination angle diameter were observed. Conversely, in the middle of the insulator, the electrical 
potential and field tend to decrease with the increase of the shed’s inclination angle diameter as can be seen 
in Figure 6. The electric field has seemed smoother with increase the shed’s radius as shown in Figure 5 and 


the inclination angle diameter as shown in Figure 6. 
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Figure 6. (a) Electric potential and (b) Electric field, with different inclination angle 


5.3. Effects of the Shed’s Diameter 

The shed’s diameter of the 11kV polymeric insulator under study is 90 mm as shown in Figure 3(c). 
The parameter is varied with values of 90 mm, 130 mm, and 170 mm to observe the effects of the electric 
field stress. The electric field stress and potential voltage along the surface of the insulator are shown in 
Figure 4. From Figure 7, it is comprehended that at the largest shed’s diameter generally yields the highest 
voltage and electric field stress [12]. However, at first shed’s outer corner, the electric field stress decreases 
as the shed’s diameter increases. 
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Figure 7. (a) Electric potential and (b) electric field, with different shed’s diameter 


5.4. Effects of the Core Diameter 

The electric potential distribution in this case has negligible changes as compared to other cases as 
illustrated in Figure 8(a). Electric field near to electrodes decreases with increasing core diameter while in the 
middle of the insulator the electric field tends to increase with increasing diameter of core insulator as shown 
in Figure 8(b). 
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Figure 8. (a) Electric potential and (b) Electric field, with different core diameter 


5.5. Effects of the Distance Between the First Shed and End Fitting 

From Figure 9 it is comprehended that the highest distance between the first shed and the bottom 
end fitting generally yields a lower electric field stress and potential voltage while the shortest distance gives 
the highest increment in electric field stress and voltage as well as. Hence, it can be concluded from Figure 
9(b) that as the distance between the first shed and the bottom end fitting increases, the electric field stress at 
first shed decreased. 
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Figure 9. (a) Electric potential and (b) Electric field, with different distance between the first shed 
and end fitting 


6. CONCLUSION 

Computation of electric field and potential distribution were carried out to study the effects of the 
various profiles on the silicone rubber insulator. The electric field stress and potential voltage along the 
polymeric insulator were highly non-uniform due to the shape of the insulator. The sharp edges of the 
insulator tend to have higher electric field stress on the insulator’s surface while the smooth edges have lower 
electric field stress. Several design considerations were studied to observe the effects to the electric field 
stress and electric potential. The design of the end fittings showed the significant change in the electric field 
stress. End fittings with round edges tend to reduce the electric field stress along the insulator’s surface as 
well as increasing the distance from sheds to electrodes. 
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